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INTRODUCTION 

Remarkable progress in recent years has resulted from the observation that the gluon 
density in a proton at small x grows as x decreases. Both DGLAP and BFKL evolution 
equations predict this rapid growth of the parton densities, thus demostrating the triumph 
of perturbative QCD. It is, however, clear that this growth cannot continue for ever, 
because it would violate the unitarity constraint [1]. Consequently, the parton evolution 
dynamics must change at some point, and a new phenomenon must come into play. As 
the gluon density increases, non-linear parton interactons are expected to become more 
and more important, resulting eventually in the slowdown of the parton density growth 
(known under the name of "saturation effect") The underlying physics can be 

described by the non-linear Balitsky-Kovchgov (BK) equation |3]. It is expected that 
these nonlinear interactions lead to an equilibrium-like system of partons with some 
definite value of the average transverse momentum kj and the corresponding saturation 
scale Qsix). This equilibrium-like system is the so called Color Glass Condensate 
(CGC) [4]. Since the saturation scale increases with decreasing of x, one may expect that 
the saturation effects will be more clear at the LHC energies. In the preasymtotic region 
{kj > Qs{x)) of energies the heavy quark production is described more adequately |0] 
by the so called ^7— factorization approach [\, 6]. At kj < Qsix) the fc^ —factorization 
approach gives us a chance to account for saturation effects (if they are under control 
the BK equation) using scale properties of the dipole model, which is equivalent the kj 
factorization f7]. 

Here we demonstrate a description of possible saturation effects in the framework of 
fcr— factorization approach. As example of that we consider the structure function at 
fixed W and low in HERA energy range and bb— quark production at LHC energies. 
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FIGURE 1. The dependence of Fl{x,Q^) (at fixed W = 276 GeV). The experimental points are 
from 1 8]. Solid curve is the result of the Arj-— factorization approach with the GLLM unintegrated gluon 
distribution from ^ . 



THEORETICAL FRAMEWORK AND NUMERICAL RESULTS 

The SF Fl at fixed W in HERA range 

The SF Fl{x, (^) is driven at small x primarily by gluons and in the fcj -factorization 
approach it is related in the following way to the unintegrated distribution 4>g(x,fc^): 

FLix^Q") = f'-f^dkl £ ej-Clix/z,Q\mlkl)^g{z,kl), (1) 

Z J i=u,d,s,c 

where ej are charge squares of active quarks. The function C|(x, Q^,mf,k^) can be re- 
garded as SF of the off-shell gluons with virtuality (hereafter we call them hard 
structure function by analogy with similar relations between cross-section and hard 
cross-section). They are described by the sum of the quark box (and crossed box) dia- 
gram contribution to the photon-gluon interaction (see, for example, Fig. 1 in To 
calculate the longitudinal SF Q^) we used the hard SF C|(x, Q^,m^,k\) obtained in 
Ref. lEo] and different unintegrated gluon distributions. ^ Notice that the fc^ -integral in 
Eq. (U]) can be divergent at lower limit, at least for some parameterizations of 4>^(jc, Z:^). 
To overcome the problem we change the low asymptotics of the QCD coupling con- 
stant within hard structure functions. We applied the so called "soft" version of "freez- 
ing" procedure flT'], which contains the shift +M^, where M is an additional 



Here we are interested in the parametrizations of unintegrated gluon functions which take into account 
sturations effects. 
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FIGURE 2. Azimutal muon-muon correlations at LHC with different unintegrated gluon distributions 



scale, which strongly modifies the infrared a, properties. For massless produced quarks, 
p-meson mass mp is usually taken as the M value, i.e. M = nipAn the case of massive 
quarks with mass m,, the M = 2m/ value is usually used. We calculate the SF as the 
sum of two types of contributions - the charm quark one F[ and the light quark one F[: 

Fl = Fl + F[. (2) 

For the F[ part we used the massless limit of hard SF and resticted ourselves to the 
modification of the argument in the strong coupling constant of the hard SF only. We 



have shown Ill2ll that our ^j-factorization results are in good agreement with the data for 
large and small part of the range. However, there was some disagreement between 
the data and theoretical predictions at 2^ ~ 3 GeV^. The disagreement comes from two 
possible reasons: additional higher-twist contributions, which are important at low 
values^, or/and NLO QCD corrections. 

It was shown that the saturation (non-linear QCD) approaches contain information 
of all orders in 1 / Q^, they resum higher-twist contributions The analysis of the 
behaviour of the longitudinal structure function Fi{x, Q^) in the saturatation models was 
done in Ref. l[l4l] ^. In Fig.l we demonstrate our A;^— factorization description of Fl{Q^) 
at fixed W with the unintegrated gluon distribution proposed in Ref. 191] which takes into 
account non-linear (saturation) effects. We see very well description of HI experimental 
data at fixed W in all region. 



Saturation effects in production at LHC 

It was shown that the data on bb azimutal correlations at Tevatron (measured as the 
decay muon ones, do /dA^^l^) are much more informative to distinguish different unin- 
tegrated gluon distributions [15]. In Fig. 2 we show our predictions for the bb azimutal 



Some part of higher-twist contributions was took into account by the "freezing" procedure. 
N.Z. thanks M. V.T. Machado for useful discussion of this problem. 
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FIGURE 3. The total transverse momentun distribution at LHC with different unintegrated gluon 
distributions 



correlations at LHC obtained with different unintegrated gluon distributions (see, for ex- 
ample, ifl^ ). The GBW and GLLM ones take into account saturation effects. We see that 
the latter two unintegrated gluon densities reduce to different behaviour of differential 
cross section Ja/JA^^^. Fig. 3 displays the total transverse momentum distribution, 
do/dpf, where pf = at LHC energy obtained with the same unintegrated 

gluon densities. In this case the difference between the results with usual unintegrated 
gluon densities and ones obtained with account of the non-linear (saturation) effects is 
more dramatic. 

CONCLUSIONS 

We considered in the framework of fc^— factorization approach possible manifestation 
of the saturation effects at HERA and LHC. We shown that account of these effects (in 
particle, by the GLLM unintegrated gluon distribution) improves to a marked degree 
description of the Fi data at low at HERA energy. We demonstrated also that the 
experimental data for bb azimutal correlations and p^^— distribution at LHC will give us 
additional possibilities to study non-linear (saturation) effects. 
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